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Emergence  of  dissipative  structures  in  current-carrying  superconducting  wires 
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We  discuss  the  emergence  of  a  spontaneous  temperature  and  critical  current  spatial  modulation  in  current- 
carrying  high-temperature  superconducting  wire.  The  modulation  of  the  critical  current  along  the  wire  on  a 
scale  of  3-10  mm  forces  a  fraction  of  the  transport  current  to  crisscross  the  resistive  interface  between  the 
superconducting  film  and  normal  metal  stabilizer  attached  to  it.  This  generates  additional  heat  that  allows  such 
a  structure  to  be  self-sustainable.  Stability  and  the  conditions  for  experimental  observation  of  this  phenomenon 
are  also  discussed. 
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I.  INTRODUCTION 

Propagation  of  electric  current  through  a  superconducting 
substance  is  an  inherently  unstable  process.  If  a  small  section 
of  the  superconductor  turns  normal,  the  heat  generated  by  the 
current  may  lead  to  a  thermal  runaway  process  (quench).  In 
practical  (stabilized)  superconducting  wires  a  normal  con¬ 
ductor  (stabilizer)  is  electrically  and  thermally  coupled  to  the 
superconducting  material.  This  allows  the  wire  to  recover 
from  accidental  temperature  perturbations  or,  if  stabilization 
fails,  to  quench  without  irreparable  damage  [1].  Thus,  a 
current-carrying  stabilized  superconducting  wire  is  a  bistable 
system  with  two  spatially  uniform  stable  modes  of  operation: 
one  is  superconducting — low  temperature  and  low  (practi¬ 
cally  zero)  dissipation,  the  other  is  normal — high  tempera¬ 
ture  (greater  than  the  critical  temperature)  and  high  dissipa¬ 
tion.  A  transition  from  superconducting  to  normal  mode 
initiated  by  a  local  temperature  perturbation  leads  to  a  nor¬ 
mal  zone  propagation  (NZP)  characterized  by  a  certain 
speed. 

Here  we  will  discuss  the  transition  from  superconducting 
to  normal  mode  of  operation  in  the  state-of-the  art 
YBa2Cu307_J.  (YBCO)  coated  conductors  [2,3].  The  main 
advantage  of  YBCO  coated  conductors  over  conventional 
low  temperature  superconducting  wires — high  operating 
temperature  (65-77  K) — has  an  undesirable  flip  side:  their 
heat  capacity  is  high  in  comparison  with  that  of  low- 7).  su¬ 
perconducting  wires.  As  the  result,  the  normal  zone  propa¬ 
gates  slowly  which  significantly  complicates  protection  of 
the  devices,  such  as  magnets  or  cables,  from  quench  induced 
damage.  We  have  investigated  a  possibility  to  increase  the 
NZP  speed  by  inserting  a  thin  resistive  layer  (interface)  be¬ 
tween  the  superconducting  film  and  the  stabilizer.  There  is 
always  a  very  thin  (a  fraction  of  a  micron)  layer  of  material 
between  the  YBCO  and  copper  that  accounts  for  the  resis¬ 
tance  to  the  current  exchange  [4,5].  Conventional  wisdom  is 
that  the  interface  resistance  has  to  be  as  low  as  possible.  If 
desired,  as  our  findings  suggest,  it  can  be  readily  increased 
by  various  means.  We  have  shown  that  the  NZP  speed  can  be 
substantially  increased  by  increasing  the  interface  resistance. 
This  is  the  result  of  the  greater  amount  of  heat  generated  in 
the  interface  during  current  transfer  from  superconductor  to 


the  stabilizer  [6].  Similar  conclusions  have  been  drawn  in  a 
large  body  of  work,  mostly  theory,  devoted  to  the  effects  of 
interfacial  resistance  on  NZP  in  conventional  superconduct¬ 
ors  [7-9]. 

The  subject  of  this  paper  is  a  finding  that  when  the  inter¬ 
face  resistance  exceeds  a  certain  threshold,  the  conventional 
scenario  of  transition  between  the  superconducting  and  nor¬ 
mal  modes  of  operation  breaks  down  and  a  current-carrying 
superconducting  wire  may  exhibit  a  much  more  complex  be¬ 
havior  than  it  is  commonly  expected.  In  addition  to  two  uni¬ 
form  states — superconducting  and  normal — there  is  also  an 
anomalous  nonuniform,  stable  mode  of  operation  character¬ 
ized  by  a  spontaneously  developing  pattern  of  spatial  tem¬ 
perature  modulation  along  the  length  of  the  wire  with  a  peak 
temperature  below  the  critical  temperature  Tc  of  the  super¬ 
conductor.  The  temperature  variation  modulates  the  critical 
current  density,  forming  the  static  temperature  and  critical 
current  density  ripples  ( T  ripples  or  Jc  ripples).  This  is  not  a 
mesoscopic,  but  a  macroscopic  phenomenon,  which  can  be 
classified  as  a  dissipative  structure — a  result  of  the  tendency 
for  the  spatially  uniform  physical  systems  driven  away  from 
thermal  equilibrium  to  break  the  translation  symmetry  and 
form  steady  macroscopic  patterns,  e.g.,  sand  ripples  [10,11], 
Taylor-Couette  flow,  thermal  convection,  oscillatory  chemi¬ 
cal  reactions,  etc.  [12].  The  spatial  scale  of  modulation  is 
determined  by  the  thermal  diffusion  length  and  in  the  case  of 
coated  conductors  it  is  on  the  order  of  0.3-1  cm,  much 
greater  than  the  magnetic-field  penetration  or  coherence 
lengths. 

II.  PLANAR  MODEL 

Coated  conductors  [2,3]  are  manufactured  in  the  form  of  a 
tape  in  which  the  superconducting  YBCO  film  of  about 
1  ptm  thick  is  deposited  on  a  buffered  flexible  metal  sub¬ 
strate  (e.g.,  Ni-W  alloy,  Hastelloy  or  stainless  steel).  A  cop¬ 
per  stabilizer  is  either  soldered  or  electroplated  on  top  of  the 
YBCO  film.  Fig.  1.  The  typical  width  of  such  a  tapelike  wire 
is  4  mm,  the  thickness,  about  evenly  divided  between  the 
substrate  and  stabilizer,  is  close  to  100  ptm  and  continuous 
pieces  with  uniform  critical  current  density  are  manufactured 
in  lengths  over  100  m  and  up  to  1  km. 
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FIG.  1.  (Color  online)  A  sketch  of  the  cross  section  of  coated 
conductor  (not  to  scale).  A  thin  superconducting  film  is  sandwiched 
between  the  copper  stabilizer  (1)  and  metal  substrate  (2).  The  resis¬ 
tive  interface  between  copper  and  YBCO  and  the  insulating  buffer 
between  YBCO  and  substrate  are  not  shown. 

In  recent  literature  the  normal  zone  propagation  in  coated 
conductors  has  been  treated  as  a  three-dimensional  (3D) 
problem  with  the  help  of  finite  element  analysis  software, 
e.g..  Ref.  [13].  In  this  section  we  show  that  the  properties  of 
the  coated  conductors  allow  the  reduction  of  this  nominally 
3D  problem  to  a  2D  or  ID  problem  without  any  significant 
loss  of  relevant  physical  content.  Computationally,  such  re¬ 
duction  makes  the  problem  much  more  tractable,  especially 
when  it  comes  to  the  analysis  of  the  effects  of  the  interfacial 
resistance  that  introduces  strong  nonlinearity  into  the  partial 
differential  equations  describing  NZR 

The  reduction  of  3D  equations  of  heat  conduction  in  a 
thin  tapelike  composite  wire  to  a  2D  (planar)  or  ID  (linear) 
model  is  possible  if  the  variation  in  temperature  across  the 
thickness  of  the  wire  is  negligible  in  comparison  with  its 
variation  along  the  wire.  The  typical  value  of  the  current  that 
coated  conductor  can  carry  is  /  — 200-300  A/cm  width.  In 
the  normal  state  the  heat  flux  per  unit  area  generated  in  the 
stabilizer 

jQ=P\J2ldu  (1) 

where  p,  and  d l  are  the  resistivity  and  thickness  of  the  cop¬ 
per  stabilizer.  At  temperatures  77-90  K  p, 

—  0.2-0. 3  puD,  cm,  and  for  <af1  =  40  pm  we  get  jg 

—  1  W/cm2.  In  order  to  remove  this  heat  through  the  surface 
of  the  conductor  a  temperature  gradient  across  the  thickness 
of  the  stabilizer  is  required: 


—  0.1-1  K,  in  comparison  with  the  variation  in  temperature 
along  the  conductor  during  NZR  This  conclusion  remains 
true  for  two  other  types  of  substrate,  Ni-W  alloy  and  stain¬ 
less  steel. 

The  last  structural  element  that  needs  to  be  examined  in 
terms  of  its  effect  on  thermal  conduction  is  the  buffer  be¬ 
tween  the  YBCO  film  and  the  substrate.  The  buffer  has  a 
total  thickness  dh  —  150-200  nm  and  consists  of  several  lay¬ 
ers  of  ceramic  substances.  If  we  take  a  representative  value 
of  the  heat  conductivity  of  yttria  stabilized  zirconia  (YSZ)  as 
Kb  —  1.5  X10-2  W/cmK,  the  temperature  drop  across  the 
buffer  that  is  required  to  transfer  the  flux  jg  is 


(4) 


Again,  as  long  as  the  thermal  flux  jg  is  on  the  order  of 
1-10  W/cm2,  the  temperature  drop  ATb  — 10-3- 10-2  K. 
The  thickness  of  the  interface  between  copper  stabilizer  and 
YBCO  is  even  smaller  and  it  is  also  can  be  considered  ther¬ 
mally  transparent.  Thus,  the  temperature  of  the  coated  con¬ 
ductors  can  be  considered  as  dependent  only  on  in-plane  co¬ 
ordinates  T(r),  r={x,y}. 

One  of  the  specifics  of  coated  conductors  is  that  the  su¬ 
perconducting  film  is  much  thinner  than  the  stabilizer  and  the 
substrate.  Therefore,  it  is  natural  to  treat  it  as  an  infinitesi¬ 
mally  thin  layer,  with  negligible  heat  capacity  capable,  how¬ 
ever,  of  carrying  finite  electric  current  density  per  unit  width 
and  generate  a  finite  amount  of  heat  per  unit  area.  In  a  tape¬ 
like  wire  shown  in  Fig.  1  we  place  the  superconducting  film 
at  z=0,  while  the  stabilizer  occupies  the  volume  0 <z<d\, 
and  the  substrate  is  located  at  -d2<Z< 0.  Since  the  esti¬ 
mates  presented  above  have  shown  that  in  coated  conductors 
the  redistribution  of  thermal  energy  and  equilibration  of  tem¬ 
perature  in  the  z  direction  is  much  faster  than  the  heat  propa¬ 
gation  along  the  {x,y}  plane,  the  two-dimensional  equation 
of  heat  conduction  can  be  obtained  straightforwardly  from 
energy  conservation.  The  energy  balance  in  a  cylindrical  vol¬ 
ume  of  height  dl+d2  and  arbitrary  cross  section  in  the  {x ,  y] 
plane  leads  to  the  following  equation: 


KXAT 

d\ 


j<2 ’ 


AT - 


Pi£ 

Kx 


(2) 


dU  ^  ? 

+  V  •  Jg  - 

dt  Q 


'  J<2,z\z=dl  j Q,^z=-d2  4  Q  ■ 


(5) 


The  thermal  conductivity  of  copper  Kx  —  5  W/cm  K,  so  that 
the  temperature  variation  across  the  thickness  of  stabilizer 
AT—  10"2-10“3  K.  A  similar  estimate  can  be  made  for  the 
substrate.  Due  to  its  relatively  high  resistance,  the  amount  of 
Joule  heat  generated  in  the  substrate  itself  is  negligible  [5]. 
The  temperature  gradient  across  the  thickness  of  the  sub¬ 
strate  is  necessary  to  transport  the  heat  generated  in  the  sta¬ 
bilizer.  Thus, 


K2AT 


AT- 


P\J2d2 

d\K\ 


(3) 


The  thermal  conductivity  of  Hastelloy  K2  —  l 
X  10“2  W/cm  K  is  substantially  smaller  than  that  of  copper. 
Nevertheless,  as  long  as  the  thermal  flux  jg  is  on  the  order  of 
1-10  W/cm2,  the  temperature  variation  across  the  substrate 
with  thickness  d2=50- 100  pL m  is  still  small,  AT 


Here  U  is  the  internal  energy  of  the  conductor  per  unit  of  its 
surface  area,  and 


where 


dU  dUdT  dT 

dt  ~  dT  dt  ~  dt  ’ 

(6) 

C  =  C\d\  +  (^2^2 

(7) 

is  the  combined  specific  heat  of  the  conductor.  The  contribu¬ 
tions  of  the  YBCO  film  as  well  as  that  of  buffer  and  the 
interface  to  the  thermal  mass  of  the  conductor  are  negligible 
in  comparison  with  that  of  the  stabilizer  and  the  substrate. 
The  in-plane  heat  flux  is  defined  as 

jQ=-KVT,  (8) 

where 
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K=K1dl+K2d2  (9) 

is  the  effective  in-plane  thermal  conductivity  of  the  conduc¬ 
tor.  The  heat  flux  from  both  top  and  bottom  surfaces  in  the 
right-hand  side  of  Eq.  (5),  we  will  take  in  the  Fourier  form 

j<2,z\z=d1  =  Kq(T-  Tq), 


stabilizer  or  substrate.  The  substrate  contributes  only  to  the 
thermal  mass  of  the  conductor  due  to  its  large  resistance  [5]. 
The  integrated  areal  density  of  heat  sources  in  Eq.  (12)  takes 
the  form 


dx  (Vi  -  V,)2  -  - 

Q=~E\+  l_s+Js-Es. 
Pi  R 


(15) 


jQ,z\z=-d2  =  K0(T-T0).  (10) 

Here  K0  is  the  heat  transfer  coefficient  across  the  insulation 
on  the  surface  of  the  wire  and  T0  is  the  ambient  temperature 
which,  in  the  absence  of  losses,  is  the  operating  temperature. 
In  the  situations  when  the  conductor  is  cooled  with  liquid 
coolant,  the  thermal  flux  from  the  surface  has  substantially 
more  complex  temperature  dependence  [4,8].  But  here  we 
are  interested  in  the  effect  of  the  interfacial  resistance,  so  it 
makes  sense  to  consider  all  other  contributing  factors  in  their 
simplest  form.  The  areal  density  of  the  internal  heat  sources 
in  Eq.  (5), 

rdi 

Q=  q(z)dz,  (11) 

J-d2 

is  the  sum  of  all  internal  heat  sources.  Thus,  the  2D  (in¬ 
plane)  heat  conduction  equations  for  the  coated  conductor 
takes  form 


Equations  (12)— (15),  supplemented  by  the  nonlinear  con¬ 
stituent  relationship  between  current  and  electric  field  in  the 
superconductor,  are  sufficient  to  solve  the  2D  NZP  problem. 


A.  ID  problem 

NZP  can  be  treated  as  a  one-dimensional  problem  when 
the  width  of  the  conductor  is  smaller  or,  at  least,  comparable 
to  the  thermal  diffusion  length,  which  we  define  below.  In 
this  case  the  thermal  equilibrium  along  the  width  of  the  con¬ 
ductor  (in  the  y  direction)  establishes  more  quickly  than 
along  the  conductor,  so  that  one  can  only  consider  the  evo¬ 
lution  of  T{x,t),  J{x,t),  and  Js{x,t )  along  the  conductor. 
Since  Jl+Js=J=  const,  the  Eq.  (13)  reduces  to  one  equation. 


dJ,  Vi  -  V, 


dx 


R 


which  can  also  be  used  in  the  form 


(16) 


-  V(KW  T)  =  Q- 2 K0(T -  T0) . 

at 


(12) 


~E  i. 


(17) 


The  redistribution  of  current  between  the  superconducting 
film  and  stabilizer  is  determined  by  the  condition  of  charge 
conservation, 

VJS  +  JZ  =  0,  VJi-Jz  =  0,  (13) 

where  Js  and  /[  are  the  linear  density  of  current  (A/cm) 
flowing  through  the  superconducting  film  and  stabilizer,  re¬ 
spectively.  The  quench  is  a  slow  process  in  comparison  to 
the  time  scale  of  charge  redistribution  in  a  metal  and  the  time 
derivative  of  the  charge  density  in  Eq.  (13)  can  be  neglected. 
The  density  of  current  flowing  across  the  interface  between 
the  stabilizer  and  superconductor 


R 

where  Vl  and  Vs  are  the  local  electric  potentials  of  the  stabi¬ 
lizer  and  superconductor,  respectively  and  A* [ £  2  cm2]  is  the 
resistance  of  the  unit  area  of  the  interface. 

There  are  three  internal  heat  sources — originating  in  the 
stabilizer,  interface,  and  superconductor,  respectively: 

q  =  —E \  + - 1 - S(z)  +  Js  ■  Es S(z)  ■ 

Pi  R 

Here  El  =  - VUj  and  Es=-Ws  are  the  electric  fields  in  the 
stabilizer  and  superconductor,  respectively.  We  use  the  delta 
functions  to  account  for  the  fact  that  two  of  the  heat  sources 
are  concentrated  in  the  volume  much  thinner  than  either  the 


The  constituent  relationship  for  a  superconductor  can  be  pre¬ 
sented  in  many  forms;  the  one  most  frequently  used  in  lit¬ 
erature  is 

Es(Jj=E0^J  .  (18) 

Here  JC{T)  is  the  critical  current  and  n  is  the  exponent,  usu¬ 
ally  large,  n~  20-40.  It  is  customary  to  take  Eq 
=  1  /xV  I  cm.  For  the  stabilizer,  the  conventional  Ohmic  rela¬ 
tionship  will  suffice: 

E\(Ji)  =  “7 -A-  (19) 

d\ 

The  one-dimensional  version  of  Eqs.  (12)  and  (15),  to¬ 
gether  with  Eq.  (16)  or  Eq.  (17),  completely  defines  the  ID 
NZP  problem.  However,  computationally,  these  two  coupled 
equations  still  present  a  fairly  formidable  problem.  Below  we 
will  introduce  an  approximation  that  allows  us  to  reduce  the 
problem  to  one  equation  that  still  preserves  the  relevant 
physics  of  the  phenomenon. 

B.  Reduction  of  ID  NZP  problem  to  one  equation 

In  good  quality  YBCO  films  the  typical  values  of  the  ex¬ 
ponent  in  Eq.  (18)  n~  20-40  [4].  We  will  define  the  internal 
heat  sources  within  the  Bean  model  approximation  [5,9] 
which  corresponds  to  the  limit  n— Correspondingly,  as 
long  as  the  current  density  Js  is  less  than  the  temperature- 
dependent  critical  current  density  Jc,  the  electric  field  and 


056224-3 

3 


LEVIN  et  al. 


PHYSICAL  REVIEW  E  79,  056224  (2009) 


FIG.  2.  (Color  online)  The  solid  line  is  a  sketch  of  the 
temperature-dependent  critical  current  density  Jc.  Temperature  in¬ 
creases  from  left  to  right.  The  dashed  line  7,  is  the  current  density 
that  flows  through  the  superconductor.  The  dotted  line  is  the  con¬ 
stant  total  current  density.  The  point  xc  demarcates  the  boundary 
between  the  subcritical  Es=  0  and  critical  Es> 0  sections. 


0=^—^’  Jc{e)=J{\-0).  (25) 

1  c~  1  1 


The  current  sharing  temperature  7]  is  defined  by  the  condi¬ 
tion  7C(7’1)=7,  so  that  7’(x|)  =  'L|  and,  therefore,  0(xi)=O. 
Equations  (23)  and  (24)  take  the  form 


A  =  J9\ 


kA  =  J 


d0 

dx 


(26) 


Thus,  the  current  in  the  stabilizer  in  the  subcritical  section 
x<xc  is  given  by 

7,  =  J8(x,)eK(x-x‘K  (27) 


and  in  the  critical  section  x>xc 


dissipation  in  the  superconductor  is  negligible.  These  we  will 
call  the  subcritical  sections  of  wire.  In  the  critical  sections  of 
the  conductor  J S=J c.  In  the  normal  state,  T>TC,  no  current 
flows  through  the  YBCO  film. 

Let  us  consider  the  situation  shown  in  Fig.  2.  The  tem¬ 
perature  rises  and,  correspondingly,  the  critical  current  den¬ 
sity  7c(x)  is  falling  from  left  to  right.  The  point  x=X\  corre¬ 
sponds  to  the  condition  J=Jc(xl).  Well  to  the  left  of  this 
point,  in  the  subcritical  section,  all  current  flows  through  the 
superconductor.  However,  in  the  vicinity  of  x1  the  current 
starts  to  leak  out  of  the  superconductor  into  the  stabilizer. 
Therefore,  the  critical  section  defined  by  the  condition  Js 
=7c(xc)  is  located  to  the  right  of  the  point  Xy  where  7 
=7f(x1).  For  x<xc  (the  subcritical  section)  Es= 0  and  the 
condition  of  charge  conservation,  equivalent  to  Eq.  (17), 


-£jy 

R—t  =  E1-Es, 

dx2 


takes  form 


-cf-Jy  pi 
R — 7  =  —J 
dx"  d\ 


1- 


The  solution  of  this  equation 


(20) 


(21) 


jl=j-jc  =  je.  (28) 

The  location  of  the  critical  point  xc  is  determined  by  the 
condition 


«e=  ±  — .  (29) 

dx 

Here  the  ±  sign  corresponds,  respectively,  to  either  positive 
or  negative  derivative  dOI  dx.  The  solution  of  the  Eq.  (29), 
0{xc)  =  0C,  is  always  positive.  Using  the  Taylor  expansion 
0(xc)~  9'(xy)(xc-xy)  we  get 

xc-xi  ~  \.  (30) 

Hereafter  a  prime  indicates  a  spatial  derivative. 

Thus,  in  the  subcritical  region  x<xc  the  heat  source,  Eq. 
(15),  is  given  by  the  sum  of  equal  contributions  from  the 
stabilizer  and  interface 

Q  =  yA  +  R(A)2  =  ^f-\(Kxc)  fe2K{x-xc\  (31) 

In  the  critical  sections  of  the  wire  where  JS=JC  all  three 
internal  heat  sources  appear — originating  in  the  stabilizer, 
interface,  and  superconductor  respectively.  From  Eq.  (16) 
follows  the  relationship  between  the  electric  fields  in  the  sta¬ 
bilizer  and  superconductor 


7,  =  AeK(x~xc\  k=  1/A  =  (pi/dyR) 1/2  (22) 

has  two  unknowns,  A  and  xc.  These  can  be  found  by  match¬ 
ing  the  in-plane  current 

A  =  J-Jc\Xc,  (23) 

and  out-of-plane  current,  see  Eq.  (13), 


J'  =  (Vi-Vs)IR,  Es  =  Ei+RJ"c.  (32) 


The  heat  source  in  the  critical  sections  x>xc  (Fig.  2)  is  given 
by 


Q=^(J-Jc)2  +  R(fc)2  +  Jc 


p± 

d\ 


U  -  7,)  +  Rf’ 


=  —  7(7-  Jc)  +  R(J’C)2  +  RJCJ ".  (33) 

dx 


dx 


(24) 


Let  us  introduce  a  dimensionless  temperature  0  and  also  as¬ 
sume  for  simplicity  a  linear  dependence  of  the  critical  cur¬ 
rent  on  temperature  [4,9] 


Expressed  in  terms  of  the  dimensionless  temperature,  Eq. 
(25),  the  heat  source  in  the  critical  section  has  the  form 

Q=—[0+\2(0')2-\2{  1-0)0"].  (34) 

dy 


Finally,  in  the  normal  sections  where  T>  Tc  all  current 
flows  through  the  stabilizer  and  the  heat  source  is  given  by 
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Q  =  —f. 
dx 


(35) 


r  > 


1 


1  -0C 


>  1, 


(41) 


In  order  to  rewrite  Eq.  (12)  in  the  standard  dimensionless 
form  and  to  reduce  the  large  number  of  material  constants  to 
their  relevant  combinations,  in  addition  to  the  definitions 
given  by  Eq.  (25),  we  will  measure  the  distance  in  units  of 
thermal  diffusion  length  lT=(DTl  y)1/2,  where  DT=KIC,  and 
time  in  units  of  y~\  where  the  increment  y=pxJ2 / dxCAT 
determines  the  characteristic  time  required  to  warm  an  ele¬ 
ment  of  the  conductor.  Here  A T=TC-TX 

Then,  piece-wise  defined  Eq.  (12)  at  T>TC  (#>1),  to¬ 
gether  with  Eq.  (35),  takes  the  form 

30  ^O  , 

—  -—2  =  \-KQ{e-0Q),  o>\, 

3t  3 


*o: 


2K0Mdx 

PxJ 2 


yt,  i=xllT. 


(36) 


In  the  critical  sections  with  temperature  0C<9<  1  Eqs. 
(12)  and  (34)  take  the  following  form: 

d  0  cP  0 

—  ~—T  =  e+  r(0')2  -  r(l  -  O)0'-KO{0-  Bo), 

3t  3(; 


0C<0<1.  (37) 

Here  tic  =  0(xc)  is  the  floating  boundary  between  the  critical 
and  subcritical  sections  determined  by  condition  (29).  The 
strength  of  the  nonlinear  terms  is  proportional  to  the  inter¬ 
face  resistance, 


r=-  =  -  R  =PA=^ZIA 
l\  Rq  °  d\  J2 


(38) 


Here  \  =  (d  tR  I  p})ul  is  the  current  transfer  length,  Eq.  (22). 

Taking  into  account  Eq.  (31)  we  get  for  the  subcritical 
regions  6<  9C, 

30  02  0  t  o  i  i 

—  -  TV 7  =  2^e-2^-Al  -  ko(0-  90),  0<  6C.  (39) 

or  oc, 

Equations  (36)-(39),  taken  together,  are  piecewise  defined 
dimensionless  version  of  the  Eq.  (12).  Although  it  is  more 
complex  than  the  standard  Kardar-Parisi-Zhang  (KPZ)  equa¬ 
tion  [14],  they  share  a  common  feature — the  competition  be¬ 
tween  diffusion  and  the  nonlinear  growth  term  °c(f?')2. 

Even  more  important  is  that  Eq.  (37)  may  have  negative 
effective  diffusion  coefficient  because  it  can  be  rewritten  as 
follows: 

d  0  0 “  0 

—  -  ( 1  -  r  +  /■  0)  — j  =  9  +  r(  O')2  -  k0(  0  -  0Q) , 


where  0C  is  determined  self-consistently  by  Eq.  (29).  Thus, 
only  when  the  interfacial  resistance,  Eq.  (38),  exceeds  a  cer¬ 
tain  threshold  we  can  expect  to  see  the  pattern  formation  and 
other  anomalous  phenomena  discussed  below. 

To  summarize,  using  the  extreme  limit  of  constituent  re¬ 
lationship  (18),  h— we  have  reduced  the  system  of  the 
two  equations,  Eqs.  (12)  and  (13),  to  one  piece-wise  defined 
Eqs.  (36)-(39)  for  temperature.  The  price  we  have  paid  for 
that  is  that  the  boundary  between  the  critical  (Et>0)  and 
subcritical  (Es=0)  regions  0,  is  a  floating  one  and  has  to  be 
determined  self-consistently  from  Eq.  (29).  This  presents  a 
significant  challenge  for  the  numerical  solution.  Our  next 
step  is  to  simplify  the  problem  by  fixing  the  floating  bound¬ 
ary  at  9C= 0.  Comparing  with  the  sketch  in  Fig.  2,  this  ap¬ 
proximation  means  that  for  x<xj  we  take  that  all  current 
flows  through  the  superconductor,  JS=J,  and  for  x>X!  we 
take  JS=JC.  This  approximation  becomes  exact  in  the  limit  of 
small  \<6lT,  see  Eq.  (30),  and  we  assume  that  even  for  X 
~  l7  the  qualitative  results  will  provide  a  meaningful  guid¬ 
ance  to  the  physics  of  the  phenomenon. 

Equation  (36)  which  describes  the  normal  section  of  the 
conductor  does  not  change.  In  Eq.  (39)  we  set  0C=O.  In  Eq. 
(37)  we  need  to  ensure  that  all  three  heat  sources  defined  by 
Eq.  (15)  are  positive.  Specifically  it  means  that  the  electric 
field  in  the  superconductor  must  have  the  same  sign  as  the 
current.  The  direction  of  the  current  flow  is  taken  to  be  posi¬ 
tive.  Therefore,  Es= pxJ/ dl(0-\20"),  see  Eq.  (32),  must  be 
positive  or  zero.  Below  we  will  use  a  step  function  to  enforce 
this  condition.  Equation  (37)  will  be  used  in  the  form 

r)  0  rP  0 

—  -—=02  +  r(0')2  +  (1  -  9){0- rff’MO- r6T) 
or  oc, 

-  K{)(0-  0O), 

0<6»<1.  (42) 

Here  the  unit  step  function  H(ES)  ensures  that  the  electric 
field  in  the  superconductor  is  either  positive  or  zero.  If  the 
boundary  between  the  critical  and  subcritical  sections  were 
determined  exactly,  this  condition  would  be  met  automati¬ 
cally.  But  since  we  set  the  boundary  fixed  at  0C=0,  the  step 
function  is  necessary  to  prevent  the  numerical  solution  to 
“wander  off’  into  unphysical  territory  where  the  electric  field 
Es  would  have  a  direction  opposite  to  current.  In  the  range  of 
temperature  that  meets  the  condition  0-r0’>  0  Eq.  (42) 
takes  the  form  of  Eq.  (37)  or  Eq.  (40).  In  the  sections  where 
H(ES)  =  0  the  heat  source  is  proportional  to  01  similar  to  that 
in  Eq.  (39). 


III.  RESULTS 


0C<0<1.  (40) 

For  r>  1  the  effective  nonlinear  diffusion  coefficient  can  be 
negative  within  the  temperature  interval  9C<  9<(r-l)l r. 
Therefore,  the  necessary  condition  for  a  negative  diffusion 
coefficient  is 


The  solutions  of  Eqs.  (36),  (39),  and  (42)  with  9C= 0  were 
obtained  by  two  different  methods  in  order  to  eliminate  the 
possibility  of  computational  artifacts.  We  used  a  direct  finite 
differences  method  and  a  high  level  software  MATHEMATICA. 
The  temperature  profiles  0(^,t)  presented  below  correspond 
to  periodic  boundary  conditions  and  the  initial  condition  in 
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FIG.  3.  (Color  online)  Sketch  of  the  right-hand  sides  of  Eqs. 
(36),  (37),  and  (39)  as  a  function  of  temperature  in  the  limit  of 
uniform  temperature.  Solid  lines  is  the  heat  source;  dashed  lines 
indicate  cooling  term  for  two  values  of  cooling  constant  k0. 


the  form  of  a  Gaussian  in  the  center  of  the  conductor 

0{t 0)  =  6>0  +  (|  0Q\  +  fl)exp{-  f/2 S2} .  (43) 

In  all  examples  shown  below  we  keep  the  width  jaf  the 
Gaussian  on  the  order  of  the  diffusion  length,  S=  \2.  The 
dimensionless  ambient  (operating)  temperature  is  determined 
by  the  ratio  of  the  transport  current  to  the  critical  current  at 
the  operating  temperature,  Eq.  (25), 

*0=l-^<0.  (44) 


The  current  sharing  temperature  is  given  by 


Tc0q-To 

0o-l 


(45) 


The  results  shown  below  correspond  to  0O=- 1  [/ 

=  ().5JC(T0)].  Correspondingly,  the  current  sharing  tempera¬ 
ture  7)  lies  halfway  between  the  operating  and  the  critical 
temperature 


Tx  = 


Tc+  Tp 
2 


(46) 


The  maximum  temperature  of  the  initial  temperature  profile, 
see  Eqs.  (25)  and  (43),  is  given  by 


Tmax=Tc  +  (a-l)(Tc-T1).  (47) 

The  results  shown  below  correspond  to  the  parameter  a 
=  1.1.  The  respective  maximum  temperature  of  the  initial 
temperature  profile  is  slightly  above  Tc.  The  stability  of  the 
conductor  with  respect  to  various  values  of  the  maximum 
temperature  was  reported  in  Ref.  [6].  The  system  described 
by  Eqs.  (36),  (39),  and  (42)  has  two  stable  uniform  ( 0'-0 ' 
=  0)  modes  of  operation,  see  Fig.  3.  One  mode  corresponds 
to  the  zero  dissipation  state  of  the  conductor  with  tempera¬ 
ture  0q.  The  other  is  a  normal  state  with  temperature 


0v -  0 o  +  *  >  1  ’  TN-T1  +  (0o  +  /q/)( Tc  -Tx). 


(48) 

The  condition  of  bistability  is 

1 

Kc\  <  Kr  —  .  . 

1  + 1 0O| 

(49) 

For  stronger  cooling,  kq  >  kc,  Eq.  (36)  does  not  have  a  stable 
uniform  solution  and  there  is  only  one  stable  uniform  mode 
of  operation — the  superconducting  state  with  temperature  0{). 

The  solutions  of  Eqs.  (36),  (39),  and  (42)  with  0C=O  fall 
into  several  categories  (scenarios)  determined  by  the  three 
dimensionless  parameters  {r,  k0,  0,}.  Here  we  present  a  few 
scenarios  that  are  important  from  both  the  basic  physics  and 
applications  point  of  view.  In  all  cases  discussed  below  we 
take  the  transport  current  J=0.5JC{T0),  which  corresponds  to 
0d=-l  and  kc=0.5.  Figure  4  shows  the  3D  plots  of  the  spa- 
tiotemporal  development  of  the  initial  Gaussian  temperature 
profile  for  different  values  of  the  cooling  constant  kq  and  the 
interface  resistance  above  the  threshold  r=  1  [negative  effec¬ 
tive  diffusivity,  Eq.  (41)]. 

In  Fig.  4(a)  a-0=0.1  (weak  cooling,  0^=9).  This  solution 
of  Eqs.  (36),  (39),  and  (42)  describes  a  conventional  bistable 
operation  of  the  stabilized  superconducting  wire.  Additional 
heat  generated  in  the  interface  (KPZ-type  growth)  substan¬ 
tially  increases  the  speed  of  normal  zone  propagation.  The 
detailed  analysis  of  the  effect  of  the  increased  interfacial  re¬ 
sistance  on  propagation  speed  and  wire  stability  is  described 
in  Ref.  [6]. 


FIG.  4.  (Color  online)  Sequence  of  temperature  profiles  as  they  develop  from  the  initial  Gaussian  perturbation  visible  in  the  foreground. 
All  profiles  correspond  to  the  interface  resistance  r=  1.4,  Eq.  (38).  The  width  of  the  initial  temperature  profile  <5=  1.4,  Eq.  (43).  The  range  of 
the  plots  are  0<yf£50  and  -40£;t//j-£40.  (a)  /c0=0.1;  conventional  normal  zone,  but  it  propagates  with  greater  speed  than  that  at  r 
<§1.  (b)  /c0  =  0.35;  conventional  normal  zone  and  temperature  ripples  emerge  simultaneously,  but  propagate  with  different  speed,  (c)  Cry- 
ostable  condition  (/c0=0.55>  kc).  Only  the  temperature  ripples  are  triggered. 
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FIG.  5.  (Color  online)  Profile  of  the  Jc  ripples  (r=1.4,  /c0=0.4) 
at  a  certain  moment  of  time,  yt= 25.  The  inset  shows  a  magnified 
view  of  the  same  profile.  The  periodicity  of  the  Jc  ripple  is  ~2.11T. 
The  peak  temperature  is  dp~ 0.24-0.26  in  agreement  with  the  exact 
value  1  -r~l,  Eq.  (51). 

At  stronger  cooling,  but  still  in  the  bistable  regime  (k0 
=  0.35  <kc),  an  anomalous  third  mode  of  conductor  opera¬ 
tion  emerges  as  shown  in  Fig.  4(b).  The  initial  Gaussian 
profile  gives  rise  to  two  distinct  fronts  propagating  with  a 
different  speed.  A  more  rapidly  propagating  front  leaves  in 
its  wake  a  metastable  state  characterized  by  the  static  spatial 
temperature  modulation  (T-ripples)  with  the  peak  tempera¬ 
ture  below  Tc.  The  temperature  modulation  causes  a  modu¬ 
lation  of  the  critical  current  density,  which  is  why  we  may 
also  call  this  pattern  7c-ripples.  The  7c-ripples  in  turn  are 
absorbed  into  the  conventional  normal  state  slowly  propagat¬ 
ing  on  top  of  them. 

When  k0>kc,  see  Eq.  (49)  and  Fig.  3,  the  cooling  power 
at  all  temperatures  is  greater  than  the  maximum  power  that 
can  be  dissipated  in  the  stabilizer  and,  as  the  result,  the  uni¬ 
form  normal  state  is  unstable.  However,  as  Fig.  4(c)  demon¬ 
strates,  for  r>  1  such  cryostable,  by  conventional  criterion, 
conductor  is  in  fact  a  bistable  system  with  the  Jc  ripples 
being  the  second  steady  mode  of  operation.  The  initial  tem¬ 
perature  profile  evolves  into  the  Jc  ripples  so  that  the  average 
temperature  of  the  conductor  behind  the  front  remains 
steady,  close  to  7)  <  Tc,  instead  of  relaxing  to  the  uniform 
superconducting  state  with  temperature  T0. 


A.  Properties  of  Jc  ripples 

In  Fig.  5  the  temperature  profile  of  the  Jc  ripples  mode  is 
shown.  At  a  distance  of  a  few  diffusion  lengths  behind  the 
propagating  front  the  temperature  modulation  is  stable  and 
alternates  in  the  range  -0.8  <  #<0.26.  In  spite  of  the  overall 
symmetry  of  the  problem,  small  computing  errors  and  inter¬ 
polation  procedures  lead  to  difference  in  the  propagation 
speed  of  the  two  fronts  and  some  irregularities  in  the  profile. 
On  the  smaller  scale,  however,  the  profile  is  rather  regular  as 
shown  in  the  inset.  The  length  of  modulation  is  2.7 lT. 

A  characteristic  feature  of  these  ripples  is  an  apparent 
discontinuity  of  the  first  derivative  SOI  dx.  This  allows  us  to 
determine  the  value  of  the  peak  temperature  0  =  max{#(x)}. 


0 


FIG.  6.  (Color  online)  The  initial  condition  gives  rise  to  two 
propagating  fronts  that  eventually  merge.  Here  r=1.5,  k0=QA ,  and 
0o  =  —  l.  There  is  a  domain  boundary  between  the  two  merged  fronts 
with  slightly  irregular  temperature  variation. 


Integrating  Eq.  (40)  over  an  infinitesimal  interval  that  in¬ 
cludes  the  temperature  peak  we  get 


(1  ~r+r0p) 


d_e_ 


=  o. 


peak 


(50) 


Since  the  derivative  is  finite,  this  condition  can  be  satisfied 
by 


0p=(r-l)lr.  (51) 

This  corresponds  to  zero  value  of  the  effective  thermal  dif- 
fusivity  in  Eq.  (40).  In  numerical  calculations  the  second 
derivative  remains  finite  and  the  value  of  9p  remains  slightly 
below  the  limit  given  by  Eq.  (51). 

The  heat  generating  regions  (#>0)  are  spaced  along  the 
wire,  segregated  by  the  colder  regions  sustained  by  heat  dif¬ 
fusion.  One  can  compare  the  average  power  dissipation  in 
the  Jc  ripples  mode  and  in  the  conventional  mode.  A  section 
of  the  conductor  between  two  minima  of  the  temperature  is 
thermally  insulated  from  the  rest  of  the  system  because  the 
thermal  flux  vanishes  at  these  points.  Integrating  Eq.  (12) 
over  the  distance  AT  between  two  troughs  we  get 

^  J  Qdx  =  2K,(Trw  -  T0),  Tav  -  ^  J  T(x)dx. 

(52) 


On  the  other  hand,  the  power  dissipation  in  the  normal  mode 

pj2 

Qn  =  =  2K0(Tn -  T0) .  (53) 

a  j 

Since  Tn>Tqv  [see  Eq.  (48)],  the  average  power  dissipated 
by  the  7,,-ripples  mode  is  smaller  than  that  dissipated  in  the 
normal  mode.  In  other  words,  in  the  regime  where  current  is 
repeatedly  exchanged  between  the  superconductor  and  stabi¬ 
lizer  the  conductor  behaves,  at  least  in  terms  of  losses,  as  a 
hyperconductor — a  hypothetical  substance  with  the  resis¬ 
tance  that  is  finite,  but  lower  than  that  of  the  stabilizer,  in¬ 
cluding  the  one  made  of  copper. 

Figure  6  illustrates  how  two  spreading  domains  of  ripples, 
originating  from  different  sources,  interact.  The  initial  con¬ 
dition  is  two  Gaussians  which  do  not  overlap.  Each  of  them 
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FIG.  7.  (Color  online)  Tristable  operation  which  corresponds  to 
/c0=0.4,  r=1.4,  and  60=- 1.  The  initial  profile  evolves  into  Jc 
ripples.  A  heater,  Eq.  (54),  is  switched  on  at  yt  =20  and  switched  off 
at  yt= 25.  It  triggers  the  transition  to  the  normal  state.  The  evolution 
of  the  system  depicted  in  the  figure  corresponds  to  £>=1.5,  8X  =  \A. 

gives  rise  to  the  spreading  T  ripples.  Once  the  spreading 
domains  merge,  a  stable  domain  boundary  characterized  by  a 
slightly  irregular  temperature  variation  is  formed. 

B.  External  heat  source  and  stability 

In  order  to  study  the  stability  of  the  ripples  pattern  we 
added  a  “heater” — an  external  heat  source  in  the  right-hand 
sides  of  Eqs.  (36),  (39),  and  (42)  that  can  be  switched  on  and 
off  at  the  moments  tx  and  t2,  respectively: 

Qex,  =  b  exp{-  e/2  8j}H(t  -  tx)H{t2  -  t) .  (54) 

Such  a  source  simulates  a  typical  experiment  on  normal  zone 
propagation  in  which  a  small  resistive  heater  is  thermally 
anchored  to  the  superconducting  wire  in  order  to  trigger  a 
transition  to  the  normal  state  [15]. 

Figure  7  illustrates  the  following  scenario.  Unlike  in  Fig. 
4(b),  the  transition  to  the  ripples  mode  is  triggered  by  the 
initial  condition  without  tripping  the  conductor  into  the  nor¬ 
mal  state.  Since  the  cooling  constant  /c0=0.4  is  below  the 
cryostability  criterion  (kq<  kc= 0.5),  the  stable  normal  state 
with  0N=  1.5  is  still  possible.  This  indicates  that  there  is  a 
finite  separation  between  the  ranges  of  stability  of  the  ripples 
mode  and  normal  mode — one  can  be  triggered  without  trig¬ 
gering  another.  At  the  moment  yt=  20  the  heater  is  turned  on 
and  turned  off  at  yt  =  25.  A  characteristic  “rhino  horn”  of 
rising  temperature  triggers  the  secondary  transition  to  the 
normal  state.  Thus,  in  this  range  of  parameters,  the  conductor 
is  a  tristable,  rather  than  a  bistable  system. 

An  ability  to  trigger  and  observe  the  temperature  ripples 
in  an  experiment  depends,  among  other  things,  on  stability 
and  robustness  of  this  pattern.  A  scenario  shown  in  Fig.  7  has 
demonstrated  a  finite  margin  of  stability  of  the  ripples  with 
respect  to  transition  to  the  normal  state.  In  Fig.  8  the  effect  of 
a  “cold  finger” — a  negative  power  source  is  shown.  Physi¬ 
cally,  this  is  equivalent  to  bringing  the  conductor  in  contact 
with  a  thermal  mass  with  a  temperature  lower  than  T0.  The 
cold  finger  was  modeled  by  a  Gaussian,  similar  to  that  given 
by  Eq.  (54),  but  offset  from  the  center  and  with  a  negative 
value  of  the  parameter  b.  In  Fig.  8  b  =  —  3,  <5j  =  8,  yf  =20,  and 


-40  -20  0  20 


x/lr 

FIG.  8.  (Color  online)  Shown  here  is  view  “from  above.”  The 
transition  to  7f-ripples  mode  is  triggered  by  the  initial  Gaussian 
profile.  The  negative  power  source  collapses  temperature  in  its  vi¬ 
cinity,  but  does  not  destroy  the  ripple  mode.  Parameters  here  are  the 
same  as  in  Fig.  6. 

yt2= 35.  After  the  application  of  the  cold  finger  the  tempera¬ 
ture  in  its  vicinity  collapses  well  below  the  background  tem¬ 
perature  0q = - 1 .  The  cold  finger  catches  and  holds  the 
nearby  propagating  front,  preventing  it  from  expanding. 
Once  the  cold  finger  is  lifted,  the  front  resumes  its  propaga¬ 
tion  with  the  same  speed  as  before.  This  demonstrates  that 
the  7- ripples  mode  is  robust  and  recovers  after  severe  dis¬ 
ruptions. 

C.  Stationary  solutions 

The  previous  scenarios  depicted  in  Figs.  4-8  show  propa¬ 
gating  fronts.  Equations  (36),  (39),  and  (42)  also  have  sta¬ 
tionary  solutions.  A  stable  temperature  profile  can  be  estab¬ 
lished  either  from  an  initial  perturbation  or  by  the  action  of 
an  external  heater.  Figure  9(a)  shows  a  simple  bimodal  struc¬ 
ture  similar  to  that  shown  in  the  inset  to  Fig.  5.  This  bimodal 
soliton  evolves  from  the  initial  Gaussian  profile.  A  more 
complex  five-peak  structure  shown  in  Fig.  9(b)  evolves  un¬ 
der  the  action  of  the  external  heater,  Eq.  (54),  with  the  same 
width  as  the  initial  condition  in  Fig.  9(a).  In  both  cases  the 
control  parameter  r=  1 .2,  Eq.  (38),  is  the  same,  as  are  the 
cooling  constant  k0= 0.5  and  the  ambient  temperature  8(l= 
-1,  Eq.  (44). 

The  stable  solitons  shown  in  Figs.  9(a)  and  9(b)  can  be 
contrasted  with  metastable,  long  lived  perturbations  shown 
in  Fig.  10.  In  these  scenarios  the  control  parameter  r  is  very 
close  to  the  threshold  r=  1.  In  Fig.  10(a)  r=  1.0426.  The  per¬ 
turbation  evolving  from  the  initial  Gaussian  profile  is  quasis¬ 
table  and  persists  for  about  yt~  22  before  dissipating.  A  very 
small  increase  in  the  interfacial  resistance  to  r=  1.0427  yields 
a  long  lived  perturbation  that  eventually  gives  rise  to  normal 
zone  propagation.  This  is  consistent  with  our  previous  obser¬ 
vation  that  an  increase  in  the  interfacial  resistance  reduces 
cryostability  of  the  current-carrying  conductor:  thermal  per¬ 
turbations  that  dissipate  at  lower  values  of  resistance  give 
rise  to  a  NZP  at  higher  resistance  [6].  What  is  important  to 


056224-8 

8 


EMERGENCE  OF  DISSIPATIVE  STRUCTURES  IN... 


PHYSICAL  REVIEW  E  79,  056224  (2009) 


7t 


50 

40 

50 

20 

10 

0 


-10  0  x/lT 


(a) 


10 


FIG.  9.  (Color  online)  View  “from  above.”  (a)  A  stationary  bi- 
modal  solution  that  evolves  from  the  initial  Gaussian  perturbation, 
(b)  A  more  complex  profile  emerging  under  the  action  of  the  heater, 
Eq.  (54).  In  both  cases  r=1.2,  /c0  =  0-5,  and  80=- 1. 


realize  comparing  Figs.  9  and  10  is  that  it  is  the  nonlinear 
KPZ  term  ~{0')2  and  the  negative  effective  diffusivity  in 
Eq.  (37)  that  allow  the  isolated  solitons  like  the  ones  in  Fig. 
9  to  be  stable  and  robust  solutions  over  a  finite  range  of 
parameters. 


D.  Thermal  oscillations 

Physical  systems  described  by  the  parabolic  partial  differ¬ 
ential  equations  do  not  support  the  propagations  of  waves.  It 
is  instructive,  however,  to  investigate  whether  the  nonlineari¬ 
ties  may  change  that  and  allow  a  long  range  propagation  of 
time -dependent  perturbations  originating  from  a  localized 
source.  The  external  power  source,  Eq.  (54)  was  modified  as 
follows: 

Qext  =  b  exp(-  £2/2<5?)[l  +  sin(2cuf)](l  -  e~'/T°).  (55) 

This  is  equivalent  to  a  resistive  heater  supplied  with  ac  cur¬ 
rent  with  frequency  a>  and  switched  on  gradually  over  a  pe- 
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FIG.  10.  (Color  online)  Shown  are  long  lived  perturbations 
evolving  from  the  initial  Gaussian  profile  with  S=  1 .4.  (a)  The  value 
of  r=  1.0426  is  very  close  to  threshold  r=l  above  which  the  effec¬ 
tive  diffusivity  can  be  negative.  The  profile  lingers  for  a  long  time, 
but  eventually  disappears,  (b)  At  slightly  greater  r=  1.0427  the  nor¬ 
mal  zone  emerges  and  starts  to  propagate.  In  both  cases  /r0=0.35, 
0O=-1,  and  <5=  1.4. 

riod  of  time  ~r0.  The  results  are  presented  in  Fig.  11. 

In  Fig.  11(a)  the  ripple  mode  is  triggered  by  an  initial 
temperature  profile  and  the  pulsating  power  source  is  gradu¬ 
ally  switched  on  with  the  time  constant  t()=  10y_1.  It  is  evi¬ 
dent  that  the  pulsations  do  not  propagate  outside  the  range  of 
localization  of  the  source.  If  there  were  a  phenomenon  of 
“thermal  sound”  it  would  manifest  itself  as  an  undulation  of 
the  peak  temperature  lines  that  run  along  the  time  line.  How¬ 
ever,  these  lines  are  perfectly  straight. 

Figure  11(b)  illustrates  how  the  confinement  of  the  ther¬ 
mal  perturbations  takes  place.  The  control  parameter  r=1.2 
and  the  cooling  constant  k0= 0.55  are  chosen  to  yield  a  bi- 
modal  stationary  soliton  similar  to  that  in  Fig.  9(a).  The  vari¬ 
able  external  power  source  creates  temperature  pulsations 
that  remain  strictly  confined  between  the  temperature  peaks. 
This  indicates  that  the  singularities  caused  by  the  negative 
effective  diffusivity  block  the  propagation  of  small  perturba¬ 
tions  through  them. 


E.  Phase  diagram 

Figure  12  presents  an  approximate  phase  diagram  that  in¬ 
dicates  the  range  of  parameters  {r,/r0}  where  different  re- 
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FIG.  1 1 .  (Color  online)  (a)  The  ripple  mode  is  triggered  by  the 
initial  Gaussian  profile.  The  pulsations  of  temperature  with  fre¬ 
quency  2 co=  1.5y  are  caused  by  the  heat  source  given  by  Eq.  (55). 
Here  r=1.5,  k0= 0.5,  60=-l,  S=Si  =  l.4,  and  b=  1.  (b)  The  initial 
Gaussian  profile  gives  rise  to  a  bimodal  stationary  soliton.  The  tem¬ 
perature  pulsations  are  confined  within  the  walls  of  the  soliton. 
Here  r=1.2,  k0= 0.55,  0O=-1,  5=  <5,  =  1.4,  and  b  =  0.3. 


gimes  take  place.  This  diagram  corresponds  to  a  particular 
value  of  00=-l  which  corresponds  to  transport  current  equal 
50%  of  the  critical  current  at  ambient  temperature  T0,  see  Eq. 
(44).  It  needs  to  be  mentioned  that  the  boundaries  between 
the  regimes  depend  not  only  on  three  parameters  {r,KQ,0o}, 
but  on  the  width  of  the  initial  perturbation  as  well.  This  is 
obvious  because  the  heat  generated  in  the  interface  is  deter¬ 
mined  by  the  derivative  of  temperature  [the  KPZ  term  in  Eq. 
(42)].  The  diagram  shown  in  Fig.  12  corresponds  to  the 
width  of  the  Gaussian  profile,  Eq.  (43),  S=  1.4. 

In  the  range  of  parameters  indicated  as  area  (1)  in  Fig.  12 
the  transition  from  superconducting  mode  to  normal  takes 
the  form  shown  in  Fig.  4(a) — a  conventional  normal  zone 
propagation,  albeit  with  increased  propagation  speed  and  re- 


FIG.  12.  The  range  of  parameters  {r,/c0}  specific  for  different 
regimes.  The  value  of  0O  =  -1  and  the  width  of  the  initial  perturba¬ 
tion  <5=1.4  are  the  same  for  all  combinations  of  {r,/c0}.  Area  (1) 
corresponds  to  conventional  NZR  In  area  (2)  T  ripples  and  NZP 
coexist.  In  area  (3)  T  ripples  are  the  second  mode  of  operation.  Area 
(4)  corresponds  to  stationary  solitons,  and  area  (5)  is  the  range  of 
cryostability. 

duced  stability  margins  in  comparison  with  that  at  low  inter¬ 
facial  resistance. 

Area  (2)  corresponds  to  simultaneous  emergence  of  two 
propagating  fronts  as  shown  in  Fig.  4(b).  Area  (3)  is  a  range 
of  relatively  strong  cooling  and  interfacial  resistance  above 
the  threshold  r=  1.  In  this  range  of  parameters  the  conductor 
is  bistable  with  T  ripples  being  the  second  stable  mode  of 
operation.  Fig.  4(c).  A  strip  (4)  corresponds  to  the  stationary 
solitons,  Fig.  9.  Area  (5)  is  the  range  of  cryostability — the 
temperature  perturbations  dissipate  without  triggering  a  tran¬ 
sition  to  the  normal  state.  Figures  10(a)  and  10(b)  illustrate 
the  evolution  of  the  initial  temperature  profile  on  both  sides 
of  the  border  between  areas  (1)  and  (5). 

IV.  EXPERIMENTAL  CONDITIONS 

In  order  to  estimate  the  interface  resistance  required  for 
the  emergence  of  Jc  ripples  we  take  the  values  of  physical 
parameters  from  Ref.  [4].  Let  us  take  the  operating  tempera¬ 
ture  7o=67  K  and  Tc~ 87  K.  At  J={M2)jf\  Tl  =  (T0 
+  TC)I2  =  11  K.  The  specific  heat  of  copper  and  substrate  at 
77  K  are,  respectively,  C^l.7  J/cm3  K  and  C2 

—  1.4  J/cm3  K.  Let  us  take  dx  =  40  /rm  and  d2= 50  pm. 
Then,  C=C\d\  +  C2d2~  14.2 X  10-3  J/cm2  K.  The  thermal 
conductivity  of  copper  and  substrate  at  77  K  are,  respec¬ 
tively,  Kl~ 5  W/cmK  and  K2~ 7  X  10-2  W/cmK.  Thus, 
K=K\d\+K2d2~203  X  10-3  W/K.  The  thermal  diffusivity 
of  the  coated  conductor  DT=KI  C~  1.4  cm2/s. 

The  self-field  critical  current  density  of  the  state  of  the  art 
coated  conductors  at  T=61  K  can  be  close  to  J{']) 

—  400  A /cm,  and,  correspondingly,  the  transport  current 
density  at  50%  capacity  can  be  7=200  Al  cm.  The  sheet 
resistance  of  the  stabilizer  is  p\ld\  —  0.5  X  10-4  ft.  Taking 
AT=10  K,  we  obtain  the  increment  y=plJ2/dlCAT 

—  14  s_1.  The  thermal  diffusion  length  lT={DTl  y)1/2 

—  3  mm.  The  characteristic  speed  of  NZP  can  be  estimated 
as 

UT=lTy~  4  cm/s.  (56) 

When  the  interfacial  resistance  is  negligible,  the  NZP 

speed  is  close  to  uT.  This  relatively  slow  propagation  speed 
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complicates  quench  detection  and  quench  protection  in  coils 
made  out  of  coated  conductors.  Increasing  interfacial  resis¬ 
tance  leads  to  increasing  NZP  speed  and  reduced  stability 
[6]. 

In  order  to  observe  the  Jc  ripples  the  interface  resistance 
has  to  exceed  the  characteristic  value  [Eq.  (38)] 

,2 

R0  =  —  «  5  X  10"6  a  cm1 2.  (57) 

d  i 

The  interface  resistance  of  the  currently  manufactured  coated 
conductors  is  about  50  nfl  cm2  [16].  Thus,  2  orders  of  mag¬ 
nitude  increase  in  interface  resistance  is  needed  in  order  to 
create  conditions  under  which  the  Jc  ripples  may  be  ob¬ 
served.  This  can  be  readily  accomplished  by  modifying  the 
stabilizer  application  procedures,  e.g..  Ref.  [17]. 

Taking  into  account  the  definition  of  the  increment  y,  the 
cooling  constant  k0  defined  by  Eq.  (36)  can  be  rewritten  as 


k0  = 


2jr()A?2/i 

Pi-/2 


2Ko 

Cy 


(58) 


To  attain  the  value  of  /c0~0.5,  see  Fig.  12,  one  needs  to  have 
the  heat  transfer  coefficient 


CyK0 

2 


«  5  X  1(T2 


W/cm2  K. 


(59) 


/;  ~  1.3  X  10“3 4 5 6 7 8 9  W/cm  K.  (60) 

Since 


h 

K0=~,  (6D 

do 

where  do  is  the  thickness  of  insulation  on  the  surface  of  the 
conductor,  it  is  sufficient  to  have  do~  200-300  /xm  to 
achieve  the  desired  range  of  k0  where  one  can  expect  to  see 
the  anomalous  phenomena  described  above. 


V.  SUMMARY 

In  summary,  we  have  uncovered  an  unusual  pattern  of 
superconducting  wire  operation  in  which  the  dc  transport 
current  is  shared  between  the  superconductor  and  stabilizer. 
This  phenomenon  manifests  itself  through  spontaneously  de¬ 
veloping  temperature  and  critical  current  modulation  along 
the  wire  with  the  spatial  scale  on  the  order  of  1  cm.  Although 
the  physics  of  this  phenomenon  is  very  different,  the  appear¬ 
ance  of  the  temperature  ripples  is  remarkably  similar  to  that 
of  the  sand  ripples.  Most  notably,  this  similarity  extends  to 
an  apparent  discontinuity  of  the  slopes  at  the  peaks  [10,11]. 
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